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It  is  demonstrated  by  a  model  computation  based  on  a  semiclasaical  theory  that  in  collisional  ionization  the 
Franck-Condon  approximation,  equivalent  to  the  assumption  that  the  ionization  event  takes  place  locally,  does 
not  always  apply. 


I.  Introduction 

We  consider  collisional  ionization  processes  of  the  type 
A*  +  B  —  A  +  B*  +  e* 

usually  referred  to  as  Penning  ionization.1  The  conceptual 
picture  most  often  used  in  the  description  of  such  processes 
hinges  on  the  Franck-Condon  (FC)  or  quasistatic  ap¬ 
proximation:  namely,  an  electron  of  energy  t  can  only  be 
emitted  at  the  particular  intemuclear  distance  ft  where 
«  -  W2(R)  -  W,(ft)  (1) 

W2(ft)  and  IV, (ft)  being  the  potential  energy  curves  cor¬ 
relating  to  A*  +  B  and  A  +  B+,  respectively.  This  picture 
of  localized  transitions,  while  intuitively  appealing  and 
generally  useful,  has  never  been  directly  verified.  Indeed, 
many  previous  approaches3'4  to  the  problem  of  Penning 
ionization  have  preferred  to  adopt  it  as  the  starting  point 
of  their  formulations.  In  this  work  we  report  on  a  model 
computation  based  on  the  He  +  Ar  system  which  dem¬ 
onstrates  explicitly  the  failure  of  the  localized-transitions 
picture.  Although  this  example  by  no  means  rules  out  the 
validity  of  the  FC  approximation  in  many  cases,10  certain 
doubts  can  now  be  raised  about  ita  indiscriminate  use. 

It  is  perhaps  useful  to  recall  the  origin  of  the  FC  ap¬ 
proximation.  The  crucial  point  is  that,  for  a  system  in¬ 
volving  bound-continuum  interactions,  the  Schrddinger 
equation  does  not  lead  to  a  differential  equation  but  to  a 
differential-integral  equation  of  the  form*  (ignoring  the 
complications  due  to  partial-wave  decomposition) 

(-  £  ^  +  "Vft)  -  E^,(£,ft)  -  /  (2) 

where 

/«  /ds/dft'  v,*(«)  v,m  P<«)  *,(£,'.«)  x 

*,(£,, vn  mem  (3) 

£0'  »  E  -  i  (4) 


(1)  For  a  raviaw  o (  thia  aubjoct,  A.  NUHaua  in  *Advanc«a  in 
Chamicai  Physic*",  Vol.  46,  Part  2,  J.  Wm.  McGowan,  Ed.,  Wilay,  Nrw 
York,  1981,  p  399  ft. 

(2)  T.  P.  O’Mallay,  Phvx.  R*v.,  If#.  M  (19##). 

(3)  M.  Mori.  J .  Phy$.  Soe.  Jpn..  24.  773  (1969) 

(4)  H.  Nakamura.  J.  Ph/$.  Soc.  Jpn.,  2#,  1473  (1989). 

(5)  W.  H.  Millar.  J.  Cham.  Phyr.  S2.  3M3  (1970). 

(6)  R.  J.  Bianiak,  Phy*.  Pao,  A,  If,  392  (1978). 

(7)  J.  C.  Bailura  and  D.  A.  Micha,  Phy*.  Pro  A,  I*.  1436  (1978). 

(8)  K.-S.  Lam,  J.  C.  Bailura,  and  T.  F.  Gaorfa,  Cham.  Pfiy*.,  U,  219 
(1978). 

(9)  J,  S.  Dahlar,  R.  R.  Tumar,  and  S.  E.  Niaiaan,  J.  Phyt.  Cham.,  84. 
1086  (1982). 

(10)  L  R..  Ralyaa  and  D  A.  Micha.  fnt  J.  (frvmtum  Cham.,  Quantum 
Cham.  Symp.,  13,  669  (1979).  In  thia  work,  tha  validity  of  tha  Franck- 
Condon  approximation  for  total  ionization  cmaa  aactiona  (aummad  over 
all  tha  partial  wavaa  and  integrated  over  •)  ia  damonatratad  indirectly 
thrmigh  tha  uaa  of  tha  two-atata  approximation  within  a  diacrattcad 
muitichannal  formaiiam.  In  our  work,  w*  attampt  to  investigate  tha 
Franck  Condon  approximation  at  tha  mor*  bmtc  level  of  ainfia  partial 
wavaa  and  ». 
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and  V,(ft)  represents  the  bound-continuum  coupling 
parametrized  by  the  electronic  energy  t  (defined  in  eq  5). 
Equation  3  clearly  indicates  the  nonlocality  of  the  poten¬ 
tial  Thia  potential,  however,  can  be  localized  with  the  help 
of  the  FC  approximation.*  It  consists  of  the  following:  in 
the  integration  over  <  in  eq  3,  only  those  values  of  c  are 
considered  important  for  the  /{'integration  at  which  the 
product  i (£#',/{ 3  ^j( EJt  0  is  slowly  varying.  If  the  wave 
functions  i(/ia  are  regarded  as  WKB  wave  functions,  and 
also  if  V.tft)  is  assumed  to  be  slowly  varying,  the  station¬ 
ary-phase  approximation,  together  with  eq  4,  would  im¬ 
mediately  lead  to  the  FC  resonance  condition  eq  1.  Thus, 
in  the  t  integral,  <  can  be  regarded  as  a  function  of  ft  (given 
by  eq  1)  in  the  product  V,*(ft)  V'.(ft').  By  virtue  of  the 
completeness  of  the  set  |^,(£,ft)|  the  rest  of  the  <  integral 
leads  to  the  delta  function  S(R  -  ft1)  and  localization  is 
finally  achieved  through  the  ft'  integration.  It  therefore 
appears  that  any  theory  which  employs  a  width  factor 
r,(ft)  *  to(»)|V,|3  in  a  complex  potential  has  already  im¬ 
plicitly  (or  explicitly)  made  use  of  the  FC  approximation: 
in  that  a  differential-integral  equation  (involving  a  non¬ 
local  potential)  has  been  reduced  to  a  differential  equation 
(involving  only  a  local,  though  complex,  potential). 

Our  approach  is  to  avoid  using  the  basis  set  of  the 
t'ectronic  Hamiltonian  (//.,)  leading  to  eq  2  and  thus  by¬ 
pass  using  the  FC  approximation  for  localization  at  this 
early  stage.  Instead,  a  basis  set  consisting  of  eigenfunctions 
of  (the  adiabatic  representation)  is  chosen,  with  the 
result  that  a  simpler  equation  consisting  of  noniocality  only 
for  «  (eq  24)  is  obtained.  This  equation  is  simpler  in  the 
sense  that  tha  kernel  £(«,«')  is  of  a  very  favorable  nature: 
it  is  separable.  We  have  shown  previously11  that  this  im¬ 
portant  feature  allows  eq  24  to  be  recast  in  the  form  of  an 
integral  equation  of  a  relatively  simple  type,  thus  per¬ 
mitting  a  solution  without  introducing  any  localization 
approximation  a  priori.  Interesting  as  it  is,  we  will,  how¬ 
ever,  not  pursue  thia  path  here.  Our  rationale  in  this  work 
ia  the  following:  having  established  a  formalism  in  which 
a  localization  approximation  ia  not  incorporated  in  the  first 
place,  we  have  the  option  of  carrying  it  through  with  or 
without  making  that  approximation  at  some  later  stage, 
and  make  a  comparison.  This  will  be  done  via  a  pertur¬ 
bative  solution  of  eq  24.  Our  purpose,  then,  is  to  inves¬ 
tigate  whether  a  localization  condition  (either  of  the  form 
of  eq  1  or  of  some  other  general  form  *  -  <(/?))  will  suffice 
for  the  solution  of  eq  24. 

Our  computation  is  based  on  a  general  semiclassical 
theory  of  collisional  ionization  developed  earlier,"  which 
treats  explicitly  the  nonlocalization  of  the  transitions  be¬ 
tween  W;  and  Wt.  In  section  II,  the  relevant  parts  of  this 
theory  will  be  briefly  reviewed  (more  complete  derivations 


(II)  K.-S.  l-.ni,  T.  F.  Gentfe,  and  D.  K.  51K.ttacH.rw,.  f'k v,  Rr<  ,-4, 

V.  13M  U9S.1). 
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Lam  and  George 
(17) 


of  equations  presented  therein  can  be  found  in  ref  11). 
Section  III  will  deal  with  the  model  computation  and  its 
results;  and  finally,  in  section  IV,  we  will  discuss  some 
pertinent  problems  for  future  investigations. 

II.  Theory 

A  diabatic  basis  set  j|02(t)),|<I'+(t),*)|  corresponding  to 
W2  and  IV,  are  assumed  such  that  the  following  conditions 


are  satisfied: 

<,fc(t)|tf|<Mt),«>  =  V,(t)  (5) 

(•MO.rWMt),*')  »  <W,(£)  +  t)&„.  (6) 

<^U)|^07«)>  *  W2(t)  (7) 

(<£>2|<A2>  =  (tb+,t\4>+,t')  =  0  (8) 

(<^2l<^+,<)  *  *  0  (9) 


with  H  being  the  total  Hamiltonian  involving  electronic 
coordinates.  Equation  5  defines  the  bound-continuum 
configuration-interaction  coupling  V,(t),  while  eq  6  implies 
that  there  are  no  couplingB  between  the  different  contin¬ 
uum  states.  (The  quantization  volume  of  the  system  will 
eventually  be  allowed  to  approach  infinity  so  that  the 
discrete  indices  will  become  continuum  ones.)  Equations 
8  and  9  specify  the  condition  of  the  absence  of  “velocity’ 
or  nonadiabatic  couplings. 

As  discussed  in  ref  1 1,  eq  5-9  are  simply  taken  to  con¬ 
stitute  the  “model’  of  the  present  theory,  and  questions 
concerning  the  explicit  construction  of  the  diabatic  states 
or  even  their  existence  will  not  be  dealt  with  here.  The 
present  theory  also  makes  use  of  classical  trajectories  for 
the  A-B  relative  motion  (with  internuclear  coordinate 
trajectory  R(t)),  so  that  V„  W2,  and  W,  all  become  explicit 
functions  of  time,  a*  they  appear  in  eq  5-7. 

We  now  choose  an  adiabatic  representation  (|d>,(f ) ) |  in 
which  nonadiabatic  couplings  do  not  vanish,  but  which 
comprise  the  exact  eigenstates  of  H: 

H|*,(f)>  -  (W,  +  <)!<£>.(£))  (10) 

These  can  be  expanded  in  terms  of  the  diabatic  set: 

|0,(t)>  =*  d,(f)|<Mt)>  +  l<^«),<>  +  Ixw(')|<Mf>.<'>  (ID 


It  can  then  readily  be  deduced  that 

-0.U)  V,.*(t) 


X„<f> 
3,(t)  =* 


-  t0+ 

-V. 


w2,u)  - ,  -  ir.u) 

with  the  instantaneous  level  width  given  by 

r,(t)  »  rpit) \Vf 


(12) 

(13) 

(14) 


where  pit)  is  the  energy  density  of  continuum  states  and 
W31  >  W,  -  W„ 

Now  |<9.(r)>  represents  Franck-Condon  “stationary" 
states,  so  that  “wave  packets’  including  the  effects  of 
nonadiabatic  transitions  between  the  |<£,(f))  can  be  con¬ 
structed: 


|^(t)>  *  Hy.it)  expj-(i/h)J^  dt'jWdf')  +  »)Jl0,(f)> 

(15) 


The  time-dependent  Sehrodinger  equation  then  yields  the 
equation  of  motion  for  y.ffh 

•y,  »  -  Lx,  (<*.k  )e  (16) 


which  implies 


7,(0)  =  /J.*(0) 


(18) 


The  transition  probability  for  the  ionization  process  is 
given  by 


where 


Jim  |S,(£)|J 


S,(f)  ^  (0+(£),<|t^(£)) 


(19) 


This  quantity  can  be  expressed  in  terms  of  the  time-de¬ 
pendent  coefficients  7,(t)  as11 

[317(0  ] 

(W21(t)  -  t)2  +  T,2(t)  J 

There  is  one  such  quantity  S,iL'  for  each  partial  wave 
characterized  by  the  orbital  angular  momentum  quantum 
number  L  of  the  A-B  relative  motion.  The  total  differ¬ 
ential  ionization  cross  section  (with  respect  to  electronic 
energy  c),  da /At,  is  then  given  by 

da/dt  -  (tP,/JV)I(2L  +  l)|S,af  (21) 

L 

where 

K0J  -  2p£,/h2  (22) 

with  p  being  the  reduced  mass  and  £t  the  total  collision 
energy  of  the  system. 

If  we  now  make  the  assumption  that  the  time  variation 
of  T,  is  insignificant,  i.e. 

If, I  «  1*2,1  (23) 


Equation  16,  on  passing  to  the  continuum  limit  in  the 
energy  variable  can  be  shown"  to  be  equivalent  to  the 
following  differential-integral  equation  for  7,: 

iT,*:, 

•  -  m 

(W21-<-,T,)1 

IVj,  V,*  r-  _  p,.v, .7, ^ 

(Wl-.  +  iTjJ  0  d<'  (Wj,  -  »'  -  T,>*  (24) 

We  note  that  the  separability  of  the  kernel  in  the  above 
equation  permits  it  to  be  recast  as  an  integral  equation — a 
Volterra  equation  of  the  second  kind."  We  will  here, 
however,  pursue  a  perturbative  solution  baaed  on  the  as¬ 
sumption 

7.  (0  -  d,-*(f)  (25) 

which  is  to  be  used  in  the  right-hand  side  of  eq  24.  Since 
3, it)  is  a  known  quantity  (eq  13),  eq  24  reduce*  to  a 
first-order  ordinary  differential  equation.  Provided  that 
for  the  collision  time*  t  of  interest 

t  «  h/2r,  (26) 


the  solution  to  thi*  differential  equation  is  gi\a?n  approx¬ 
imately  by" 


1  r 1 

7,U)  *  7.(9)  +  - —  J  df ' - - 

2irp,)'nJ°  r,1 


,Ci/*)(.-W„>lV«r./* 


5(»  -  Wj,  -  ,T.l 


The  integral  is  of  a  form  such  that  the  stationary- phase 
approximation  ia  likely  to  be  applicable  for  its  evaluation. 
This  approximation  would  be  exact  in  the  strict  classical 
limit  5  *0;  and  in  a  great  variety  of  collision  problems 


The  initial  condition  ir  given  by 


Co— lonal  lonbatton  u  a  Nontocaizsd  Process 
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Figure  1.  PotartW  cuvoa  tor  *w  Ha  +  Ar  co-ton*  tonfcslton  system. 
W+P)  correlates  to  He'lltSs,^)  +  Ar  aiv.  WJfl)  to  Ha  +  Ar".  Al 
quantities  an  in  atomic  units 


signs  on  the  incoming  and  outgoing  portions  of  the  tra¬ 
jectory,  the  “true*  stationary  configurations  (/?< r0)>  may 
be  shifted  to  either  side  of  the  Franck-Condon  configu¬ 
rations  by  an  appreciable  amount. 

At  this  point,  we  have  no  a  priori  reason  to  expect  that 
the  stationary-phase  approximation  will  work  for  eq  27. 
Indeed,  it  may  be  surmised  that,  since  the  present  problem 
involves  bound-continuum  interactions,  localization, 
whether  it  conforms  to  the  Franck-Condon  points  or  not, 
would  not  be  expected  to  take  place.  In  the  ner  sectic  n 
we  will  attempt  to  settle  the  question  by  perfon  ling  cal¬ 
culations  on  a  specific  model  system. 

III.  Model  System  and  Results 

We  apply  eq  27  to  the  He*(ls2s,3S)  1-  Ar  Penning  ion¬ 
ization  system  using  the  following  model  potential  curves 


(all  values  are  givsn  in  atomic  units):8 

WJR)  -  Vj (R)  Ola) 

W,(fl)  *  VX(R)  -  E0  £0  »  0.149  (Sib) 

V,(R)  -  A,*'*1*'  -  8t{R)  *,(*)  »  •  1,  2  (32) 

8i(R)  •  Cgi/R*  +  Cn/R*  (33) 


MR)  -  [1  +  expHR-RJAini  +  Ml  +  exp((R  - 

Rod/bi m  (34) 


treated  in  the  semiclassical  context  which  do  not  involve 
bound-continuum  interactions,  it  is  indeed  found  to  be 
excellent.11  The  advantage  of  the  stationary-phase  ap¬ 
proximation,  when  valid,  ia  that  the  result  allows  a  very 
intereating  physical  interpretation  of  the  collision  process, 
namely,  transitions  take  place  locally  in  configuration 
space.  In  the  context  of  eq  27,  the  lowest  order  station¬ 
ary-phase  result  would  be  given  by 


7  .**<*>  =« 

7.(0)  + 


1 _ Y' 

,-UTM j/s  (  ;  )  I 

V°"  "W/ 


(28) 


In  this  equation  the  sum  is  over  the  stationary  points  f0 
which  satisfy 

«  -  W„(t)  -  Wn(t)t  (29) 


the  ±  signs  in  the  exponential  are  to  be  used  according  sa 
-2Wh  -  W„i  >  or  <  0.  If  eq  28  ia  valid,  the  most  im¬ 
portant  contributions  to  the  ionization  procaaa  would  then 
be  localized  to  the  discrete  times  to,  or  the  discrete  con¬ 
figurations  R(tg),  where  R  is  the  A-B  intemuclear  coor¬ 
dinate.  The  fact  that  there  may  be  more  than  one  sta¬ 
tionary  point  for  a  particular  trajectory  leads  to  quantum 
interference  effects  within  a  semiclassical  treatment 
We  recall  that  in  most  conventional  interpretations  of 
Penning  ionization,  localization,  according  to  the  Franck- 
Condon  approximation  (eq  1),  is  assumed  to  take  place 
where 

i  -  W„(t)  -  0  (30) 

Equation  29,  however,  indicates  that,  when  localization 
applies  at  all,  it  takea  place  at  points  quite  different  from 
those  implied  by  eq  30.  In  fact,  since  Wtl  has  opposite 


(12)  W.  H.  Millar,  Adu.  Chtm.  Phyt..  21.  #9  (1974). 


A,  -  44.65270 

A,  -  4.39678 

(35) 

6,  -  0-63.10 

62  »  0.9675 

(36) 

C„  -  7.94 

C„  -  226.0 

(37) 

C4l  -  0.6904 

C«  -  0.0 

(38) 

-  6.42317 

0,  -  2.81505 

(39) 

Rm  -  5.76899 

R„  -  8.73010 

(40) 

These  potential  curves  are  shown  in  Figure  1  (drawn  to 
scale).  The  interaction  configuration  coupling  strength  T, 
is  taken  to  have  a  constant  value: 

I\  -  1.066  x  10“*  (41) 

This  assumption  of  constant  coupling  strength  may  be 
somewhat  unrealistic  as  T,  is  usually  Uken*  to  have  an 
exponential  dependence  on  R.  It  is  adopted  here  to  render 
the  model  as  simple  as  possible,  in  particular,  as  the  va¬ 
lidity  of  localization  does  not  depend  on  the  nature  of  T,. 

Computations  are  made  for  d^/dr  (zeroth  partial  wave) 
for  a  range  of  values  of  t  between  0.166  and  0.19  au  ac¬ 
cording  to  eq  21  and  27.  For  a  particular  value  of  «,  the 
classical  trajectory  is  obtained  by  propagation  on  the  po¬ 
tential  curve  W,  +  e,  using  Hamilton's  equations  of  motion. 
The  choice  of  the  potential  surface  follows  as  a  conse¬ 
quence  of  eq  10.  The  initial  kinetic  energy  E,  for  the 
trajectory  corresponding  to  a  given  <  ia  then  given  by 

E,  *•  Et+  E0  -  t  (42) 

where  Ex  is  the  actual  initial  kinetic  energy  with  respect 
to  the  potential  curve  Wj(f?).  For  the  present  computa¬ 
tions  this  is  chosen  to  be 

Ex  -  0.061  (43) 

The  reduced  mass  of  the  He-Ar  system  used  here  is  given 
by 

M  -  6631.406  (44) 

Propagation  ia  initiated  at  the  intemuclear  distance  R  » 
20.0  in  the  approeching  phase  and  carried  out  until  R  — 
20.0  in  the  separating  phase.  Results  are  obtained  for  (i) 
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Figure  2.  WJt  and  IVj,  a  a  functions  of  time  for  total  Initial  kinetic 
energy  F.  =  0  »ui  and  t  =  0.19.  71  quantities  are  In  atomic  units. 
Results  are  obtained  py  propagating  Ha  mat  on's  equations  wtoi  a  time 
step  size  St  =  10  »■;. 

the  direct  numerical  trapezoidal  quadrature  of  eq  27,  in 
which  case  the  time  step  size  for  the  propagation  is  de¬ 
creased  until  the  integrated  values  become  convergent,  and 
(ii)  the  stationary-phase  case  using  eq  28,  in  which  case 
the  time  step  size  is  taken  to  be  10  au.  (The  totai  time 
required  for  the  above-specified  trajectory  is  ~  1.3  X  10s.) 
Figure  2  shows  the  results  for  W21  and  Vv2,  as  a  function 
of  time  for  <  =0.19  au  and  propagation  stepwise  SC  -  10 
au. 

In  Figure  3,  the  results  are  presented  for  d<rl0,/d<  (<  = 
0.155  to  0.19  au)  using  the  two  methods  of 
computation — direct  numerical  quadrature  and  station¬ 
ary-phase  approzimation.  These  results  demonstrate 
clearly  the  discrepancy  between  the  two  methods,  indi¬ 
cating  the  failure  of  the  stationary-phase  approximation, 
and  consequently  invalidating  the  picture  of  localized  in¬ 
teractions.  It  is,  however,  interesting  to  note  the  quali¬ 
tative  resemblance  between  the  two  seta  of  results. 

IV.  Discussion 

The  main  conclusion  to  be  drawn  from  this  work  is  that 
the  Franck-Condon  approximation  (in  the  sense  stated  in 
the  introduction)  cannot  be  applied  indiscriminately  to 
bound-continuum  problems,  such  as  Penning  ionization. 
This  implies  that  the  convenient  picture  of  localized 
transitions  in  configuration  space  does  not  necessarily  hold 
in  molecular  collisions  involving  a  continuum  of  energy 
eigenstates.  We  should  note,  however,  that,  even  though 
strict  localization  (at  a  discrete  configuration)  has  been 
invalidated  “rough"  localization  may  still  hold  because  of 
the  factor  (<  -  Wn  -  iT,)'1  occurring  in  the  integrand  of 
eq  27.  But  then  the  degree  of  localization  is  much  harder 
to  specify  than  in  the  stationary-phase  case.  The  strongest 
statement  one  can  make  is  that,  since  i  -  W2]  -  if,  -  0 


Figure  3.  L  =  0  partial  wave  differential  Ionization  cross  section 
dih°Vd <  as  a  function  of  <.  The  dashed  curve  is  the  stationary-phase 
result,  obtained  by  using  eq  25;  the  solid  curve  gives  the  exact  result, 
obtained  by  numerical  quadratire  of  the  integral  in  eq  2-4. 

roughly  corresponds  to  the  Franck-Condon  condition  eq 
1,  the  ionizing  transition  rate  peaks  at  the  Franck-Condon 
points.  The  extent  of  the  actual  “spread"  of  the  transition 
region  remains  quite  indefinite.  Hence,  the  Franck-Con¬ 
don  approximation,  though  still  useful  as  a  crude  estimate 
of  the  general  vicinity  of  the  transition  region,  cannot  be 
used  as  a  basic  premise  justifying  strict  localization  in 
collision  theory.  Moreover,  in  cases  where  the  station¬ 
ary-phase  approximation  does  apply,  strict  localization 
obtains,  and,  as  shown  in  section  II,  the  Franck-Condon 
approximation  will  predict  erroneous  localization  points. 

Our  mode!  for  the  theory  has  been  chosen  to  illustrate, 
under  the  most  simplifying  assumptions  possible,  the  in¬ 
adequacy  of  the  Franck-Condon  approximation.  As  a 
consequence,  some  compromise  to  reality  has  been  made. 
One  such  is  the  assumption  of  the  constancy  of  T,  with 
respect  to  R.  Since  the  validity  (or  lack  of  it)  of  the 
Franck-Condon  approximation  should  not  depend  on  the 
detailed  nature  of  T,,  this  assumption  has  been  made  in 
the  interest  of  arriving  at  the  relatively  simple  form  for 
the  differential-integral  equation  of  motion  (eq  24)  and 
the  subsequent  result  eq  27.  These  would  require  modi¬ 
fication  if  the  time  (or  spatial)  variation  of  f,  becomes 
important.  Hence,  further  work  including  this  effect 
should  prove  worthwhile  in  order  to  bring  the  model  closer 
to  reality,  although  the  major  conclusion  is  not  expected 
to  be  changed. 

Another  worthwhile  direction  of  further  investigations 
would  be  the  nonperturbative  treatment  of  eq  24  via  the 
solution  of  an  integral  equation  (as  mentioned  in  section 
II).  or  the  carrying  out  of  the  perturbative  solution  to 
higher  orders.  It  must  be  pointed  out  that  we  have  not 
attempted  here  to  delineate  rigorously  the  general  criteria 
of  validity  for  the  perturbative  treatment  used.  This 
problem  would  also  appear  to  warrant  further  work. 

With  regard  to  the  physical  characteristics  of  the 
problem,  a  natural  extension  is  to  include  the  spin  and 
propagation  direction  of  the  emitted  electrons  as  observ¬ 
ables.  This  would  entail  enlargement  of  the  continuum 
basis  set  to  include  these  degrees  of  freedom,  and  also  some 
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detailed  studies  of  the  spin-  and  direction-dependent 
configuration-interaction  coupling  strengths.  Finally,  it 
is  interesting  to  point  out  that,  when  these  extra  degrees 
of  freedom  are  included,  the  problem  becomes  formally 
identical  with  that  of  molecular  fluorescence  in  a  collision 
context,13  in  which  the  photon  frequency  <■>,  wave  vector 
direction  k,  and  polarization  direction  <  correspond  to  <, 
k  (electronic  momentum  unit  vector),  and  S  (electronic 

(13)  K.-3.  Um  and  T.  F.  Gtorgt.  J.  Chtm.  Phys.,  7*,  3396  (1982). 


spin  unit  vector)  in  the  problem  of  collisional  ionization. 
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